Nonlinear evolution of the plasma heatwave: Compressing the laser beatnotes via 

electromagnetic cascading 
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The near-resonant heatwave excitation of an electron plasma wave (EPW) can be employed for 
generating the trains of few-femtosecond electromagnetic (EM) pulses in rarefied plasmas. The 
EPW produces a co-moving index grating that induces a laser phase modulation at the difference 
frequency. The bandwidth of the phase-modulated laser is proportional to the product of the plasma 
length, laser wavelength, and amplitude of the electron density perturbation. The laser spectrum is 
composed of a cascade of red and blue sidebands shifted by integer multiples of the beat frequency. 
When the beat frequency is lower than the electron plasma frequency, the red-shifted spectral 
components are advanced in time with respect to the blue-shifted ones near the center of each laser 
beatnote. The group velocity dispersion of plasma compresses so chirped beatnotes to a few-laser- 
cycle duration thus creating a train of sharp EM spikes with the beat periodicity. Depending on the 
plasma and laser parameters, chirping and compression can be implemented either concurrently in 
the same, or sequentially in different plasmas. Evolution of the laser heatwave end electron density 
perturbations is described in time and one spatial dimension in a weakly relativistic approximation. 
Using the compression effect, we demonstrate that the relativistic bi-stability regime of the EPW 
excitation [G. Shvets, Phys. Rev. Lett. 93, 195004 (2004)] can be achieved with the initially 
sub-threshold heatwave pulse. 

PACS numbers: 52.35 Mw, 52.38 Kd, 52.38 Bv 



I. INTRODUCTION 



An electron plasma wave (EPW) is a natural tool for 
manipulating the properties of intense radiation beams. 
It can be used for up-shifting the laser frequency m , for 
enhancing the self-focusing of co-propagating 0, Q and 
counter-propagating 4] radiation beams^for the resonant 
self- modulation of the laser amplitude 0, and for cou- 
pling the signal and pump lasers in the parametric ampli- 
fier [6(. Also, the high-amplitude EPW driven by a short 
laser pulse can induce the pulse shrinkage with time 7]. 

Excitation of the EPW by the ponderomotive force 
(beatwave) of the two-color laser with the difference 
frequency Q close to the electron Langmiur frequency 
ujp = y/ 4ire 2 no/m e (no is an electron plasma density, 
m e and — |e| are the electron rest mass and c harg e) has 
attracted attention for a long time || |j| 0, The 
nonlinear plasma wave is highly sensitive to the varia- 
tions of frequency and amplitude of the ponderomotive 
force. For example, by chirping the beat frequency [p3 ] 
the EPW excitation can be enhanced by the autoreso- 
nance effect. Downshifting the beat frequency from the 
plasma resonance (f2 < ui p ) can also result in the large- 
amplitude wake excitation due to the effect of relativistic 
bi-stability 0, 0] . Conversely, the driven electron den- 
sity perturbations can cause the laser amplitude modula- 
tion, either transverse 0, Q or longitudinal 0] . There- 
fore, the performance of the beatwave scheme critically 
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depends on the self-consistent evolution of the light and 
plasma waves that includes effects of numerous plasma 
nonlinearities 0|. The relativistic self-phase modula- 
tion i stimulated forward Raman scattering 
(SFRS), and electromagnetic cascading I5L ITrL 12(1 EHj 
(EMC) broaden the laser frequency bandwidth, while the 
group velocity dispersion (GVD) of radiation distorts the 
laser amplitude. In certain regimes this scenario results 
in a strong local enhancement of the laser field. Specifi- 
cally, the frequency downshifted (Q < u> p ) beatwave pulse 
of initially low amplitude can be transformed into a train 
of sharp electromagnetic (EM) spikes of relativistic in- 
tensity and few-laser-cycles duration [l^ . The spikes are 
separated in time by the beat period tj, = 27r/f2. 

In the present paper a theoretical model is formu- 
lated which, in a weakly relativistic approximation, accu- 
rately describes the nonlinear excitation and relativistic 
bi-stability of the EPW, and frequency and amplitude 
modulation of the laser pulse in one spatial dimension 
(ID) and in time. We select specific regimes in which 
the EMC induced by near-resonantly driven EPW causes 
laser spectral broadening. At every point of the per- 
turbed plasma, the EPW creates an index grating co- 
moving with the laser beams. Hence, the periodic fre- 
quency modulation (FM) of the laser develops at a dif- 
ference frequency O. In spectral terms, the FM is man- 
ifested as a cascade of Stokes and anti-Stokes sidebands 
shifted by integer multiples of f2 from the laser funda- 
mental Wo- 

The effect of EMC has been known in physics of laser- 
plasma interactions since early 70s when Cohen et al. poj 
suggested to enhance plasma heating by the decay of the 
cascade-driven EPW. Later on, the EMC was considered 
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FIG. 1: Schematic of the two-stage cascade compressor. 
The frequency modulation (FM) occurs in a rarefied plasma. 
Denser plasma is used for the compression of beatnotes. 



as an EPW diagnostic in the plasma heatwave accel- 
erator [2l|. Systematic study 0, Q of plasma waveg- 
uiding options provided by the nonlinear interaction of 
laser beams with the cascade-driven EPW revealed an 
enhanced self-focusing of the co-propagating beams de- 
tuned in frequency below plasma resonance (f2 < w p ). 
Calculations of Refs. 0, describe the non-stationary 
cascade evolution in two dimensions (2D) in the planar 
and three dimensions (3D) in cylindrical geometry, take a 
full account of the relativistic nonlinearities of both cas- 
cade components and the EPW, but neglect the GVD 
of radiation and, hence, the longitudinal (temporal) dis- 
tortion of the laser amplitude caused by the longitudinal 
transport of the cascade energy. Our work fills this gap 
by concentrating on the ID compression of laser beat- 
notes due to the EMC and GVD. 

We prove that the GVD completely dominates the evo- 
lution of weakly-relativistic heatwave (with intensity of 
initial beams over 10 16 W/cm 2 ) in either centimeter-scale 
rarefied (no ~ 10 18 cm -3 ) or millimeter-scale dense plas- 
mas (no > 10 19 cm -3 ). We show that the plasma wave 
driven below the plasma resonance, < lu p , chirps the 
laser frequency in a very special way: near the center 
of each laser beatnote the red-shifted sidebands are ad- 
vanced in time with respect to the blue-shifted ones. The 
GVD can compress thus chirped beatnotes to a few-laser- 
cycle duration provided the laser bandwidth tends to ojq. 
The effect of GVD is controllable: proper adjustment 
of the plasma and laser parameters can reduce it while 
preserving the desirable bandwidth. In this case, the 
cascade compression can be made in two stages [H (see 
Fig. Pi : (i) a low-density plasma (the Modulator) with 
fl < LUp(M) is used for the FM of initially two-frequency 
laser, and (ii) a higher-density plasma (the Compres- 
sor) with Wp(c) 3> ^ serves for the beatnote compres- 
sion. Therefore, a train of sharp electromagnetic spikes 
of intensity by several orders of magnitude higher than 
ionization threshold for any medium can be generated. 
Similar concept of using Raman cascades for radiation 
beams compression in molecular gases has been success- 
fully tested in experiments |22j at low laser intensities 
(<C 10 14 W/cm 2 ). The technique of the EMC compres- 



sion in gases is not appropriate for applications requiring 
high laser intensity. One such application is using multi- 
ple short laser pulses with a tunable time delay for the co- 
herent generation of plasma waves [2^, l24l |25| . Making a 
sequence of several independent ultrashort high-intensity 
laser pulses with the periodicity of less than one picosec- 
ond could be a major experimental challenge 23]. The 
approach discussed in our paper suggests a viable path to 
creating such pulse trains at weakly relativistic intensity. 

The outline of the paper is as follows. In Section ITTlwc 
derive the basic theoretical model (subsection 111 Aj l and 
analyze the basic scalings for the EMC and the cascade 
compression (subsection III Bll . In a realistic plasma, the 
EMC is a complicated interplay between the sideband 
coupling through the driven electron density perturba- 
tions, GVD of radiation, nonlinearities due to the rel- 
ativistic increase of an electron mass, and SFRS. Fully 
nonlinear simulations presented in Sec. 1 1 1 1 1 account for 
all these effects and describe the cascade development in 
either two-stage (subsection IIII All or single-stage (sub- 
section lTlI B|) compressor. Because the longest time scale 
of the problem is only a few ion plasma periods, para- 
metric decay of the EPW [2(| and consequent plasma 
heating |20j are insignificant and thus ignored. The sim- 
ulation parameters of subsections IIII Al and IIII 51 are op- 
timized so as to make the relativistic nonlinearities and 
SFRS almost negligible. When the parameters of the 
setup are not optimized, the SFRS can be seeded by the 
plasma wake driven by the heatwave pulse of finite du- 
ration. Contribution from the SFRS into the cascading 
process is discussed in the subsection llll CI In the subsec- 
tion IIII Dl we show how the beatnote compression helps 
in the EPW excitation via relativistic bi-stability [jjl 
with the initially sub-threshold laser intensity. Conclu- 
sion gives the summary of the results. In Appendix lAlthe 
amplitude of the plasma wake excited by a given detuned 
heatwave pulse is evaluated. 



II. ONE-DIMENSIONAL THEORY OF EMC 



A. Basic equations 



We assume that the laser duration does not exceed a 
few ion plasma periods, so the ions are immobile and form 
a positive neutralizing background. In one spatial dimen- 
sion and in the limit of weakly relativistic electron mo- 
tion, Maxwell's equations and hydrodynamic equations 
of electron fluid give the coupled equations for the longi- 
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tudinal and transverse momentum of electrons [2j 
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Here, a = p eJ _/(m e c) and q z = p ez /(m e c) are the nor- 
malized transverse and longitudinal components of the 
electron momentum, and q 2 = q 2 + a^ < 1. We take 
a = Re (eoa), where eo = (e x + ie y )/y/2 is a unit vector 
of circular polarization; hence, q 2 = q 2 + \a\ 2 /2. In the ID 
approximation, conservation of the transverse canonical 
momentum expresses the normalized momentum through 
the laser vector potential, Aj_ = (m e c 2 /e)a. The normal- 
ized electron density perturbation, 
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obtained through differentiating Eq. Hla(l with respect to 
z and t. At the plasma entrance 2 = 0, the amplitude of 
a planar two-frequency laser beam is given by 

a(0,t) = e-^ nt [ao(0,<) + ai (0,<) e ^], (4) 

where Q rs uj p -C Wo- The ponderomotive beatwave [the 
right-hand side (RHS) of Eq. J3J] produces an electron 
density grating co-moving with the laser beams. The 
moving index grating produces the cascade of laser side- 
bands, 
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where ka — Cl/v g , and v g — k$c 2 /u>q is the group ve- 
locity associated with the laser fundamental frequency. 
v g is found from d = n /n c = 1 — {v g /c) 2 , where 
n c = m e ujQ / (Aire 2 ) is a critical plasma density. The am- 
plitudes a n vary slowly in time and space on scales fi -1 
and fc^ 1 . 

To describe the nonlinear evolution of the cascade J^J , 
we take into account Eq. and rewrite Eq. (|lbl) as 
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We retain in Eq. JBJ) the terms of order not higher than a 3 . 
Having in mind that q z ~ Sn/riQ, and that the relativistic 
saturation of the beatwave-driven EPW occurs at q z ~ 
a 2 / 3 , we keep in Eq. © the nonlinear terms of order a 3 , 
ack, aq 2 , ag 3 and finally arrive at 
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originating from the relativistic mass correction of an 
electron oscillating in the transverse (R a ) and longi- 
tudinal (R q ) fields describe the relativistic self-phase- 
modulation of laser. The leading nonlinear current term 



C = a(Sn/no) 



(9) 



is responsible for the EMC and stimulated forward Ra- 
man cascade |28|. Our earlier work 0] assumed the 
non-resonant EPW excitation (Sn/no ~ a 2 ), so the 
term Il8b(l was neglected. By including this term we in- 
clude the regime with relativistic saturation of the reso- 
nantly driven EPW (Sn/n ~ a 2 / 3 ^> a 2 ). 

We substitute the expansion (J5J into Eq. (JJJ, replace 
the variables (z,t) by (where £/v g = t — z/v g is a 

retarded time, and z is the propagation distance through 
plasma), and collect the equal frequency terms. The re- 
sulting set of coupled envelope equations 
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a n « k 2 (C n -R a n -R n ), (10) 



where k p = lo p /c, accounts for the propagation of side- 
bands through plasma [the first term in the left-hand side 
(LHS)], the GVD of the sidebands (the second term in 
the LHS), and the sideband coupling through the non- 
linearities (the RHS). To evaluate the RHS of Eq. fHty 
we have to specify the nonlinear plasma response. The 
nonlinear electron density perturbation is driven by the 
ponderomotive force [the RHS of Eq. J3J|] approximated 
as -{c/2) 2 Y J i{lkn) 2 Pi{z^y lkni - Here, 
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(ii) 



and \dpi/dz\ <C fco|p/|. We expand the density perturba- 
tion in the ponderomotive force harmonics, 



(12) 



where <5rt_; = 5n*, \5ni\ <C rio, and \d6ni/d£\ -c fefj | <5rji | , 
and assume that each density harmonic is driven by the 
corresponding harmonic of the ponderomotive force. In 
the expansion (|12f> . the terms with I = ±1 are the closest 
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in frequency to the natural modes of plasma oscillations. 
They produce the dominant contribution to the cascade 
dynamics. Keeping in Eq. © the terms of order not 
higher than a 4 , and having in mind the scaling |<5n/no| ~ 
a 2 / 3 that holds in the case of relativistic saturation of the 
resonantly driven EPW Q, we find that the amplitude 
N e = 6ri—i(z,£)/n c obeys the nonlinear equation 



(13) 



Here, 5ui = (Q 2 — Wp)/(2fi) is the beatwave detuning 
from the plasma resonance, and R is proportional to the 
nonlinear frequency shift due to the relativistic mass in- 
crease of an electron oscillating in the longitudinal and 
transverse electric fields, 
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The initial condition for Eq. 1)1 3|l is N e (z, — oo) = (un- 
perturbed plasma ahead of the pulse). Amplitudes of 
the non-resonantly driven density harmonics (I =/= ±1) 
are found from the linearized Eq. @: 
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Using equations l|15fl . we evaluate the nonlinear terms 
in the cascade equations I|1U[) . We extract from the terms 
C a contribution from the EPW harmonics of orders I ^ 
±1 and include it into the terms R a . Therefore, only the 
contribution from the near-resonant EPW harmonic iV e 
determines the form of the "cascade" nonlinearity C n . 
And, only near-resonant EPW harmonic is taken into 
account for evaluating the terms R^. The result is 
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The second term in the RHS of Eq. (|16b(l comes from 
the nonresonant EPW harmonics (|15J) . prime means that 
the terms with I = ±1 are not included in the sum. The 
physical meaning of the nonlinearities l|16|) is as follows. 



• The terms C n couple the neighboring laser side- 
bands through the near-resonantly-driven har- 
monic of the EPW. These terms describe the elec- 
tromagnetic cascading and the stimulated forward 
Raman cascade. 

• The terms i?^ describe the nonlinear frequency 
shifts produced by the relativistic mass increase of 
electron oscillating in the transverse fields and by 
the nonresonantly- driven harmonics of EPW. 

• The terms Rl describe the nonlinear frequency 
shifts produced by the relativistic mass increase 
of electron oscillating in the longitudinal electric 
field of the near-resonantly-driven harmonic of the 
EPW. Only when the EPW is driven resonantly 
and reaches the relativistic saturation the term R^ 
can dominate i?^. In all the simulations that will 
follow in this paper i?^'s are negligibly small. 

Assuming v g « c, we rewrite the set (1101) as 
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The boundary condition for Eqs. i|17fl is given by Eq. 10}. 
Equations l|13|l and l|17|) form the basis of ID spatio- 
temporal weakly relativistic theory of the near-resonant 
plasma beatwave excitation. This theoretical model en- 
compasses the phenomena of the nonlinear excitation and 
relativistic bi-stability of the EPW 0[I3, periodic FM 
of the lase r Il5l 1221. GVD of the cascade components, 
and SFRS [iMillll. 

By the judicious choice of parameters the terms C n 
can be made dominating in the RHS of Eq. (TJJ, and 
laser spectral broadening will occur exclusively due to 
the EMC. Despite the large bandwidth of laser achieved 
in certain regimes of EMC, the effect of GVD can be neg- 
ligible (see the discussion at the end of subsection III Bll . 
The laser amplitude can be then modified in a separate, 
denser, plasma with the high GVD (the Compressor). In 
this two-stage scenario, nonlinearities of the Modulator 
affect primarily the laser phase, while in the Compressor 
the GVD modulates the amplitude. Compression of the 
laser beatnotes in plasma can result in the laser intensity 
so high as to give |et| 2 ~ 1. 

The laser sidebands in Compressor remain coupled 
through the nonlinear frequency shifts, and the laser fre- 
quency bandwidth keeps growing. Thus, the relativis- 
tic nonlinearities of plasma can compete with the linear 
compression process. The Compressor density is cho- 
sen so as to entirely exclude the possibility of resonant 
plasma response: u p (c) is never close to an integer mul- 
tiple of fl. Provided \np n \ <§; 1, the amplitude of den- 
sity perturbation at the n th beatwave harmonic is de- 
scribed by Eq. (|15l) with no and u> p replaced by n (c) 
and ^ p (c)- As the electron density perturbations in the 
Compressor are nonresonant and thus small, we neglect 
the terms i?^. Moreover, the terms C„ are absorbed by 
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R® [that is, summation in the second term of Eq. l|16bl) 
is extended to I = ±1, and C„'s do not show up in the 
Compressor equations] . We redefine the retarded time as 
C/ v g(C) = t- z / v a (C) (where v g{c) is the group velocity 
of the laser fundamental component in the Compressor 
plasma of density roo(C) 3> no) and find that the com- 
pression process can be described in terms of the coupled 
nonlinear equations similar to Eqs. (|17|) : 



2i d , 

Kq OZ 



u 



^0 

0J„ 



-dc^K^. (18) 



Here, dc = no(c)/ n c <C 1 is the normalized Compressor 
density. The boundary conditions are given by the solu- 
tion of Eqs. I|17|) at the Modulator exit, a n (z = zm>0- 




FIG. 2: Normalized distance between the points of the first 
maximum and half-maximum of Jm(z) (solid line) and the 
scaling function M~ 2 ^ 3 (dashed line). 



B. Basic scalings for laser frequency modulation 
and compression 

In the ideal two-stage compressor, the processes of 
EMC and compression are separated. The EMC devel- 
ops in the Modulator plasma with zero GVD, while in the 
dense Compressor plasma with all the nonlinearities ne- 
glected the GVD compresses the radiation beatnotes. It 
is instructive to derive the basic scalings for each process 
because these approximate scalings will help to select the 
optimal parameters of fully nonlinear simulations. 

When both R^ and R^ are taken to be zero, and the 
GVD is neglected (d — 0) in Eq. I|17|) . scaling laws for 
the EMC are particularly simple. Assuming u> n ~ wo, 
we derive from Eqs. H17|) a set of conservation laws: 
dpi/dz = 0. Hence, pi = for / ^ 0,-1, Po{ z ,0 = 
M(U)| 2 + M0,£)| 2 , p-i(z,t) = ai(O,0o8(O,O. and, 
in the co-moving frame, N e is independent of z despite 
the evolution of the laser phase. Thus simplified Eqs. I|17|) 
have the analytic solution [2lJ 

a»(*,0= E ^(O^y^^Jn^W), (19) 

CT=0,1 

satisfying the initial condition [here, J n {x) are the 
Bessel functions, and ijj(z,£) and W(z, £) are the phase 
and absolute value of the generating function w(z,£) = 
We^ = i(k z/A)N e (^)}. Substituting Eq. JT^J into 
Eq. JSJ yields the expression for a train of phase- 
modulated beatnotes, 

a (z,0 = zl a n (0,£)cos[k n t + ip(z,£)}, 

rc=0,l 

where tp(z,£) = [k z /2)\N e (£)\ sin(V> - fc n f). Tne physi- 
cal meaning of this result is that, without GVD, the laser 
undergoes frequency modulation only. The magnitude of 
the plasma wave depends only on the laser amplitude 
which remains unchanged. This is valid for any pair of 
ao(0, £),ai(0, £) and the corresponding N e (£). 

The FM is periodic in time with the beat period 
Tb when iV e (£) is almost constant [this is the case for 



\dp-i/d£\ <C \(5u/c)p-i\]. To avoid the oscillations of 
N e with time due to the relativistic dephasing 0, 0, 0] , 
we take \Su\ > 3(w p /4) \/3|p_i| 2 /2. Then, the term pro- 
portional to 5ui dominates in the LHS of Eq. |JT3}, which 
yields N e (£) « d[p_i(£)/4](Q/<5u;); then, for real p-i(£), 

a= ]T a„(0, cos [fc„£+ (k z/2)N e (0 cos(fc^)] • (20) 

n=0,l 

From equation l|19|l . a Modulator plasma slab of thickness 

z M » 2A4/(V e fc ) (21) 

produces .M sidebands on either side of the fundamental, 
and a frequency bandwidth Aw ~ 2^/dMujQ. Conversely, 
M. r~j r e ZM^o\n e — n-o\, where Ao = 2-itc/ujq is the fun- 
damental laser wavelength, and r e = e 2 /(m e c 2 ) is the 
classical electron radius. 

As follows from Eq. I|20|) . only for 5uo < the laser 
wavelength decreases with time near the amplitude max- 
imum of each beatnote (positive chirp). The GVD of 
plasma tends to compress thus chirped beatnotes: the 
shorter (blue-shifted) wavelengths catch up with the 
longer (red-shifted) wavelengths, eventually building up 
the field amplitude near the center of each beatnote. 
Thus, a sequence of sharp spikes is produced. If we con- 
sider an unperturbed Compressor plasma of given den- 
sity, neglect the relativistic effects by setting R^i' w 0, 
and fix laser frequency and the number 2M. of satellites, 
we find that the peak compression occurs at a distance 

zc * (^) (£) ■ (22) 

This estimate assumes that the outer sidebands were ini- 
tially separated in time by roughly Tf,/2 within one beat- 
note. To catch up with the red sidebands at the beatnote 
center, the blue sidebands need the propagation time 
zc/c w (c/ Av g )(Tb/2), where the group velocity mis- 
match is Av g w 2Mil(dv g /dui) UJo w (3Mfl/k )(ujp/LLiQ) 2 . 
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FIG. 3: The two-stage cascade compression. The physical 
quantities are shown at the entrance (z = 0, light gray) and 
exit of the Modulator [z = Zs, medium gray), and after the 
Compressor (z — zg + zc, black), (a) The laser pulse intensity 
(the time window contains about 100 beatnotes). (b) The 
normalized amplitude of the near-resonant EPW, Sn-i/no = 
N e /d. (c) The beatnote intensity near the laser pulse center; 
one beat period near £ = is shown, (d) The laser spectra 
near £ = 0. The nonlinearities and GVD in both plasmas are 
included. 

The nonzero GVD of radiation in the Modulator 
plasma must be properly accounted for. The cascade 
components can be redistributed in time and space thus 
reducing coherence of the EPW excitation and affect- 
ing the frequency chirp. Naively, the GVD can be- 
come significant in the Modulator whose length zm is 
close to the compression length estimated from Eq. <|22[) 
with w p ((7) = io p . However, the higher-order Stokes-anti- 
Stokes sidebands are generated later in plasma and have 
less time to catch up with the fundamental. Recalling 
that \cim\ oc |Jm(z)| in the Modulator, we define the half- 
growth length z M i/ 2 at which \clm\ reaches a half of its 
maximum value, | Jm( z m,i/2)\ — \Jm{ z m)\/^- Thereby, 
compression effectively takes place over the shorter dis- 
tance Az w zm — z M .i/2 < zm- The analytic formula 
Az w M~ 2 / 3 zm accurately fits Az for M > 5 (see 
Fig. Therefore, Az <C z M for M » 1. Hence, for 

z M < M 2/3 z c pa (\ /6)M- 1/3 d- 2 , (23) 

the effect of the GVD is negligible in the Modulator be- 
cause the distance Az actually available for the compres- 
sion is less than zc- Otherwise, if zm ^ Ai 2 ^ 3 zc, the 
GVD in the Modulator becomes important. 

Another manifestation of the GVD in the Modulator is 
the SFRS seeded due to finite duration of the beatwave 



pulse. The stimulated forward Raman cascading |28( can 
interfere with the process of phase modulation and con- 
taminate the laser frequency chirp. Reduction in the 
compression efficiency can follow. The effect of SFRS 
is examined in subsection IIII Gl 

III. NONLINEAR SIMULATIONS OF THE EMC 

A. The two-stage cascade compressor 

We model the EMC by numerically solving the set 
of coupled nonlinear equations (|13(l and (|17|l with the 
boundary condition 

a (0,£) = ai((U)=^-« 2 A c ^ 2 (24) 

for the laser sidebands, and N e (z, £ = — oo) = for the 
EPW. The beatnote compression in the second stage is 
modeled by numerically solving the Compressor equa- 
tions H18|) . All the nonlinearities associated with the 
effects of rclativistic mass correction and non-resonant 
electron density perturbations are retained in the mod- 
eling of both stages. In all the simulations below, the 
fundamental laser wavelength is fixed at Ao = 0.8 [im. 

The two-stage compression starts with the initial laser 
amplitude A — 0.2, the Modulator density no = 8.75 x 
10 17 cm~ 3 (hence, d = 5 x 10~ 4 ), and 5u) = -0.1w p . 
The laser pulse duration is tl = 4.5 ps (about half 
the ion plasma period for a fully ionized Helium). Hav- 
ing chosen the maximum density perturbation |iV e (z = 
0)|max ~ 0.5 x 10~ 4 and resulting spectral width of the 
laser (M « 8 sidebands on each side), we find the Mod- 
ulator length z% ~ 4.1 cm (such interaction length could 
be implemented in a plasma channel |29j). 

The simulation results are shown in Fig. |3J From the 
plot (a) it is seen that the peak laser intensity at the 
Compressor exit (z = zg + zq) is by a factor of 7.2 larger 
than at the Modulator entrance [z = 0). The increase 
in intensity results from the shown in plot (c) beatnote 
compression from the initial duration of T^ in ^ ~ 120 fs 
to r h ( ottt -) w 13 fs (roughly 5 laser cycles). Compres- 
sor plasma has the density n (c) = 25no and length 
« 0.0275zg sa 1.1 mm (such a short dense plasma can 
be created by ablation of a microcapillary |30l|). 

The inequality l|23|) is very well satisfied for the Mod- 
ulator parameters. Consequently, the beatnote pre- 
compression seen in the plot (c) is quite insignificant. 
Plot (b) shows that N e {z,£ t ) also reveals almost negligible 
variation with z in the Modulator. Thereby, according 
to the plots (b) and (c), the EMC develops in accordance 
with the scenario outlined in subsection III Bl 

Compression in the second stage clearly proceeds in 
the nonlinear regime. The laser amplitude becomes rel- 
ativistic (\a\ — * 1), and the nonlinear frequency shifts in 
Eqs. I|18l) couple the laser sidebands and further increase 
the laser bandwidth. Figure Old) shows that the result- 
ing frequency spectrum is at least twice as broad if com- 
pared with that at the Modulator exit. As a consequence, 
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the linear formula l)22[l overestimates zc by a factor of 
three since it ignores both pre-compression of the pulse 
in the Modulator and additional bandwidth increase in 
the Compressor. Also, quality of the compressed beat- 
notes is not perfect: instead of a single sharp spike, one 
can observe a multi-spike structure in Fig.^c), the dis- 
tance between the spikes being roughly T&/5. One can 
relate this structure to the phase modulation occurring 
due to the electron density perturbation at fifth harmonic 
of the heatwave frequency, which is the closest to the nat- 
ural mode of the Compressor plasma oscillations. Due to 
this effect, one beatnote is not gradually compressed into 
one spike but rather splits into five spikes, of which the 
one located near the original beatnote maximum has the 
largest amplitude. 

A number of issues are yet to be addressed before the 
theory of the cascade compression is complete. The ne- 
glected effects of transverse evolution of the laser, such 
as relativistic self- focusing and cascade focusing [3j, are 
dominant at high plasma density in the Compressor and 
are potentially adverse. But, we find the ID scenario 
of the two-stage compression conceptually simple and 
useful for understanding the underlying phenomena. In 
the next subsection we consider the single-stage approach 
which assumes concurrent cascading and compression in 
the same low-density plasma. 



B. The single-stage cascade compressor 

Increasing plasma density in the Modulator increases 
the GVD. Therefore, we can explore the idea of com- 
pressing the beatnotes concurrently with generating the 
sidebands. In the following set of simulations the elec- 
tron density is doubled, no = 1.75 x 10 18 cm~ 3 . Cascade 
compression is simulated in plasma of the same length, 
z s = 4.1 cm, and with the same initial density pertur- 
bation \N e (z — 0)| max s» 0.5 x 10~ 4 as in the previ- 
ous subsection. The laser amplitude is A w 0.071, and 
the heatwave detuning is Soj = — 0.025w p . As we shall 
see in subsection IIII CI the SFRS manifestation can be 
large in this regime. Appropriately low seed level can be 
achieved with the heatwave pulse envelope varying slowly 
on the time scale Sui^ 1 . We choose the heatwave pulse 
tl = 14.25 ps long (which is about three ion plasma pe- 
riods for a fully ionized Helium), which corresponds to 
|5u;tl| ~ 27. The given initial intensity on axis and laser 
duration require the pulse energy of about 5 J in the focal 
spot of 30 /im radius. 

The linear estimate of the effective compression 
length (122[) shows that the inequality (|23|l almost breaks 
under the simulation parameters. Hence, a beatnote 
compression is large at the plasma exit. Figure Ufa) 
shows that the resulting peak intensity is by an order 
of magnitude larger than at the plasma entrance. Ini- 
tial duration of the laser beatnote, as shown in Fig.^fc), 
is reduced from r b ( in ) w 85 fs by roughly a factor of 10 
(to roughly 3 laser cycles). The spectral features of the 
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FIG. 4: The single-stage compressor with concurrent EMC 
and compression. In plots (a) - (c) , the physical quantities are 
shown at the entrance [z = 0, gray) and exit of the plasma 
(z = z$, black), (a) The laser pulse intensity (the time win- 
dow contains about 500 beatnotes). (b) The normalized am- 
plitude of the near-resonant EPW, <5n_i/no = N e /d. (c) The 
beatnote intensity near the laser pulse center, £ = 0. (d) The 
laser spectra near £ = with (black) and without (gray) GVD 
and all nonlinear frequency shifts. 



EMC are different from those obtained with the nonlin- 
ear frequency shifts and the GVD neglected [i.e., with 
d = in Eqs. i|17fl ]. A red asymmetry of the cascade 
spectrum is seen in Fig.^fd). Some spectral broadening 
versus the case of d = can be attributed to the self- 
phase-modulation produced by the nonlinear frequency 
shifts. Figure Bib) shows that the electron density per- 
turbation is not an integral of motion. Its final amplitude 
is roughly twice the initial, and a small-amplitude wake 
is left behind the train of compressed spikes which can be 
recognized as a signature of SFRS. Despite of this, nei- 
ther the self-phase-modulation nor the SFRS are adverse 
for the cascade compression. 

The transverse evolution of the cascade is a matter 
of high importance for the experimental verification of 
this compression scheme. For instance, electron density 
perturbations can significantly lower the nonlinear fo- 
cusing threshold of counter-propagating laser beams Q . 
The co-propagating cascade of electromagnetic beams 
also experiences enhanced focusing in both 2D (planar) 
and 3D (cylindrical) geometry if SI < ui p (Refs. 0, 0)- 
We shall give here a few necessary estimates (effects of 
transverse evolution will be given a detailed consider- 
ation in the upcoming publications). When the beat- 
wave is downshifted by Suj — — 3(w p /4) ^/3|p_i| 2 /2 the 
self-focusing threshold in the planar 2D geometry 0, 
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is ao(fcpXo/10) 3 > 0.064, where xq is a laser focal spot 
size (ao = aoe~ x l x °). So, the sub-threshold regime un- 
der the parameters of Figs. [3] and 0] requires the spot 
size Xq < 40 /im. If we loosely translate Xq into the 
radius of the laser focal spot in the cylindrical geome- 
try, we find that the refraction-limited interaction length 
2zr = (2tt/X )xq is less than 1.2 cm. However, the re- 
quired zg 4.1 cm can be achieved by means of the 
plasma channel guiding j2(| . 

In the regime of cascade compression considered above, 
eliminating potentially adverse effects of laser self-phase- 
modulation and SFRS required complying with some 
hard restrictions on the laser pulse amplitude, dura- 
tion, and beatwave frequency detuning (\a\ 2 <C 1, and 
|(5a;rL| ^s> 1). The next set of simulations shows that 
these conditions are desirable but not necessary for the 
manifestation of the effect. The cascade compression 
can be observed even for ultrashort (~ 100 fs) beatwave 
pulses propagating in a dense plasma (no ~ 10 19 cm -3 ) 
where neither GVD nor relativistic nonlinearities are 
small. We consider the evolution of a two-color ultrashort 
laser [3ll | whose energy is initially distributed between 
the fundamental (97%) and the Stokes (3%) components. 
The laser frequencies are u>o = 2.356 x 10 15 s _1 (Ao = 
0.8 fjxa) and lo x = 2.159 x 10 15 s" 1 (A x = 0.873 urn). As- 
suming that loq — u>\ = 0.95o; p , we derive the plasma den- 
sity n = 1.35 x 10 19 cm" 3 ; hence, d = 7.725 x 10~ 3 . We 
choose ao « 0.3, and ai « 0.048. At z = 0, the laser pulse 
is Gaussian l)24|) with a duration tx « 90 fs. In this case 
|(Jwrt| ~ 0.9, and the beatwave pulse amplitude is not 
slowly varying. Nevertheless, under these seemingly un- 
favorable conditions, the EMC develops very effectively 
and the intensity contrast of the amplitude-modulated 
laser pulse grows rapidly. The plasma length is chosen 
so as to produce 5 sidebands on either side; the plasma 
length is bounded from above by z§ m 2.5 mm, while 
the compression length evaluated from formula i)22|) for 
M. = 5 gives the lower bound, z rs 0.45 mm. The most 
spectacular features of the EMC shown in Fig. [3] are ob- 
served atzK 1.8 mm. At the plasma border, z = 0, the 
intensity variation of the two-color laser is about 50%, 
while at z « 1.8 mm very deep amplitude modulation 
develops with the intensity contrast ratio reaching a fac- 
tor of 25. The mostly affected are the beatnotes in the 
tail of the laser pulse; they are compressed to roughly a 
quarter of a plasma period. The laser nonlinear evolu- 
tion boosts the amplitude of the plasma wake (which is 
increased by a factor of 3). The laser spectrum broadens 
and reveals a red shift by about — u>o/2. This is a clear 
indication of the forward stimulated Raman cascade [2^ 
which, as appears in this simulation, does not prevent 
the beatnote compression. Moreover, the laser frequency 
red-shifts towards the pulse tail. The red-shifted field 
components that form the compressed beatnotes in the 
tail move slower than those in the pulse head. Thus, as 
seen in Fig. EJa), the beatnotes in the tail accumulate 
a considerable time delay (about a quarter of the beat 
period) with respect to the initial positions of their max- 
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FIG. 5: The ECM of the two-color short-pulse (90 fs) 
laser [3l| in a dense (1.35 x 10 19 cm -3 ) plasma. Physical 
quantities are shown at the plasma entrance (z = 0, gray 
lines) and exit (z = 1.8 mm, black lines). 



ima. This time delay is in agreement with the frequency 
shift —uio/2- 

The simulations presented in this subsection demon- 
strate the robustness of the EMC in the conditions when 
the GVD is large, and the nonlinear processes of the rela- 
tivistic self-phase modulation and the SFRS interfere the 
cascading process. Slow variation of the beatwave pulse 
envelope is therefore helpful but not necessary for the 
cascade compression. 
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FIG. 6: Electron density perturbation (a) and the temporal 
profile of laser intensity (b) for the parameters the same as 
of Fig. ^except the laser duration reduced by a factor of 2.5, 
Upi~L ~ 100. Gray color — z — 0, black — z = z%. 



FIG. 7: Electron density perturbation (a) and the temporal 
profile of laser intensity (b) for the parameters the same as of 
Fig- El except the laser duration reduced to lu p tl « 440, and 
|<5cjtl| = 10. Gray color — z — 0, black — z = z$. 



C. Manifestation of SFRS in cascade compressor 

Equations (|17|) admit the longitudinal transfer of elec- 
tromagnetic energy in the co-moving frame. Therefore, 
the laser amplitude modulation may result not only from 
the EMC with concurrent compression of beatnotes but 
also as a consequence of the SFRS instability 0, 0, |2£| 
(also referred to as the ID resonant modulational insta- 
bility 5\). The SFRS is different from the EMC. The 
latter is merely a phase modulation which may proceed 
in the absence of GVD. 

The SFRS is a resonant process seeded by the electron 
density perturbations oscillating at the plasma frequency 
Dp. The instability bandwidth is much narrower than oj p 
even in the case of relativistically strong pump, ao ~ 1 
(Ref. |l9|). And, in the examples of Figs. and 01 the 
SFRS bandwidth is much lower than the absolute value 
of the beatwave frequency detuning Su>. However, an 
electron plasma response to the laser beatwave always 
contains a component oscillating at u> p . This component 
is due to the finite duration of the beatwave pulse, and 
its amplitude is governed by the product |Jwrx,|. These 
resonant density perturbations can be enhanced by the 
SFRS to a level comparable to that of a non-resonant 
plasma response, and can interfere the phase modulation 
process. Hence, the effect of the SFRS is adverse and 
should be avoided by the judicious choice of laser and 
plasma parameters. 

The seed level for the SFRS can be reduced by taking 
|(Jwtl| 3S> 1. For example, parameters used in Figs. [31 
and 0] correspond to |<Wtl| ~ 30 and reveal no SFRS 
manifestation: no considerable plasma wake is left behind 
the laser at z = zg. Hence, the plasma response is almost 
entirely non-resonant in these simulations. 

Under parameters of Fig. reduction in the beatwave 
pulse duration by a factor of 2.5 (|<Wti,| = 10) pro- 
duces visible enhancement of the plasma wake that can 
be attributed to the SFRS manifestation (see Fig. |HJ). 
At the plasma border, the wake amplitude is 5n(z — 
0) = Sn s fa 2.4 x I0- 3 n . Taking Sn s as the SFRS 
seed amplitude, we can theoretically evaluate the am- 



plification factor by using the formula (4.12) of Ref. [l9j . 

\n\5n(z)/5n s \ ps (2/c) (z /+~ r§(r) dr) 1/2 . This expres- 
sion takes into account the laser temporal profile and 
is valid for z > £ > cr L ; rg( T ) = (l/8)(w*/a;g)a§(T) 
stands for the instantaneous growth rate. The theoretical 
estimate of amplification factor is In \6n(zg) / 5n s \ ps 3.36. 
On the other hand, comparing the wake amplitudes at 
the entrance (z = 0) and at the exit (z — z s ) of the 
plasma gives the amplification factor of In \8n(zg) / 5n s \ ps 
3.2, which is very close to the analytical estimate. We 
have found that the theory and simulation agree for 
80 < ujpTL < 300. Therefore, throughout this range, 
plasma wakes are excited almost entirely by the SFRS. 
Remarkably, the maximum laser intensity, as well as the 
shape of individual beatnotes, is almost the same at 
z = zg, for the parameters of Figs. |3 and Hence, in 
the considered parameter range the effect of SFRS has a 
negligible effect on the laser evolution. 

Oppositely to the just discussed case of rarefied 
plasma, reducing the laser duration by the same factor 
2.5 under the parameters of Fig. 0] (i.e., plasma twice as 
dense versus that of Figs. and |HJ) causes significant en- 
hancement of SFRS. Figure Efa) shows the plasma wake 
amplification by a factor of In \8n(zs)/8n s \ « 6.05. Theo- 
retical estimate of the SFRS gain gives 4.2; this discrep- 
ancy can be partially explained by the laser amplitude 
growth due to the beatnote compression. Figure 0Tb) 
demonstrates the strong deformation of the beatwave in- 
tensity profile. The spectral content of the electromag- 
netic cascade varies considerably and can exhibit cither 
overall red- or blue-shift at different £ in the window 
-100 < k n £ < 100. 

In conclusion, to get rid of the SFRS, one should keep 
the SFRS seed low by keeping the product |<5cjtx| large. 
As the simulations show, it should be larger than 20; 
this may require the beatwave pulse duration of several 
picoseconds or larger. 
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D. Relativistic bi-stability of EPW 

Cascade compression is a perfect tool for studying 
threshold phenomena. One of them, the relativistic bi- 
stability (RB) of the EPW driven by the long (\5u>t l \ > 
1) beatwave pulse with fl < u> p , is considered in this 
subsection. The RB results in the excitation of large- 
amplitude plasma wakes 0|. The intensity threshold 
should be met for the RB to occur. The threshold is 
multi-faceted: it is determined by the beatwave frequency 
detuning 5oj, the laser amplitude, shape, duration, contri- 
bution from the plasma wave harmonics etc. Various as- 
pects of the RB in the approximation of the given driver 
P-i{C) are addressed in the forthcoming publication [l4| . 

Evolution of the laser amplitude becomes important 
when a laser propagates over a considerable distance in 
the plasma. Strong distortion of the beatwave tempo- 
ral profile within a finite distance (few millimeters) in a 
dense (no ~ 10 19 cm~ 3 ) plasma does have adverse conse- 
quences for the amplitude- and phase-sensitive process of 
RB. On the other hand, we demonstrate below that the 
beatnote compression helps to cross the RB threshold in 
the case of initially sub-threshold laser amplitude. We 
start the simulation with the parameters of the numer- 
ical example from Rcf. 13] : simulation starts at z = 
in a plasma with a density uq = 10 19 cm~ 3 , the beat- 
wave pulse having a Gaussian temporal profile Ij24(l with 
ujpTL ~ 212 and 5u> — — 0.05w p , and the laser fundamen- 
tal wavelength being Ao = 0.8 fim (then, d « 0.0057). 
Solution of Eq. i|13|) with a given driver (equivalent to 
the calculations of plasma response at the entrance point 
2 = 0) show that the RB threshold is A 2 m 0.03375 (16% 
lower than in Ref. IT3lh This threshold corresponds to 
the normalized peak intensity |a|^ ax « 0.135. To demon- 
strate how this threshold is crossed in the course of laser 
evolution in plasma, we start the simulation with a sub- 
threshold value of the laser intensity, A 2 w 0.027. The 
simulation results are shown in Fig. [SJ 

As the laser travels through plasma, the cascade com- 
pression of the beatnotes locally increases the inten- 
sity, and the RB threshold is crossed at a distance 
z w 0.53 mm. At this point, the magnitude of the elec- 
tron density perturbation jumps abruptly by a factor of 
2.3 [Sri-i/no immediately before and after crossing the 
RB threshold is shown in Fig.|SJa)]. After that point, the 
resonant density perturbation grows steadily to roughly 
80% of background density. At z w 1.2 mm the beat- 
wave amplitude distortion becomes so strong as to de- 
stroy the coherence of the plasma response, and <5n_i 
drops sharply. Importantly, at the point where the RB 
threshold is met, the beatnote amplitude is not very dif- 
ferent from sinusoidal [see the inlay in Fig. 18Tb)]. and 
the shape of the beatwave pulse is not much different 
from the initial Gaussian. The normalized intensity at 
which the RB occurs in the simulation is |a|^ ax ~ 0.145. 
The amplitude of the beatwave pulse immediately be- 
fore (z = 0.52 mm) and after (z = 0.54 mm) crossing 
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FIG. 8: Relativistic bi-stability of the EPW. Magnitude of 
the near-resonant electron density perturbation [plot (a)] and 
the laser intensity [plot (b)] are shown at z — (light gray), 
z = 0.52 mm (medium gray), and z = 0.54 mm (black). Inset 
in plot (b): one beatnote selected near £ = 0. Plot (c): the 
laser spectrum at z = 0.54 mm. The RB threshold is crossed 
at 2 w 0.53 mm. Crossing the threshold increases the EPW 
amplitude by a factor of 2.3. 



the RB threshold is almost the same, as can be seen in 
the plot IHtb). As follows from Fig. [Sfc), the beatnote 
compression necessary for reaching the RB threshold is 
achieved at the laser bandwidth roughly equal to Wo/3. 

Therefore, the proposed model of EMC in plasmas with 
nonzero GVD is able to demonstrate the effect of rela- 
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tivistic bi-stability in the dynamic simulations with ini- 
tially sub-threshold laser amplitude. 

IV. CONCLUSION 

In this paper, we have developed a nonlinear model 
that describes the evolution of laser beatwave and elec- 
tron density perturbations in time and in ID in space in 
the weakly relativistic regime. Electromagnetic spectrum 
evolution and the effects of finite group velocity disper- 
sion are accurately modeled. The model includes the 
nonlinear frequency shifts related to the relativistic cor- 
rections of electron mass and the harmonics of the elec- 
tron density perturbations. It also takes into account the 
spatio-temporal evolution of the near-resonantly driven 
electron density perturbation. The theoretical model also 
describes a number of nonlinear effects important for the 
implementation of plasma beatwave accelerator. It is 
found that, for the beatwave downshifted in frequency 
from the plasma resonance, the electromagnetic cascad- 
ing produced by the near-resonant electron density per- 
turbations leads to the compression of the laser beat- 
notes, which finally transforms the beatwave pulse into 
a train of sharp (few-laser-cycle) electromagnetic spikes 
separated by the beat period in time and space. A train 
of electromagnetic pulses useful for the particle accelera- 
tion applications [23, 13 E3 can be self-consistently cre- 
ated. We are also able to demonstrate how the electron 
plasma wave of large amplitude can be excited due to the 
effect of relativistic bi-stability even in the case of initially 
sub-threshold beatwave pulse. The work is supported 
by the U.S. Department of Energy under Contracts No. 
DE-FG02-04ER54763 and DE-FG02-04ER41321, by the 
National Science Foundation grant PHY-0114336 admin- 
istered by the FOCUS Center at the University of Michi- 
gan, Ann Arbor. 

APPENDIX A: ELECTRON DENSITY 
PERTURBATION DRIVEN BY A GIVEN 
BEATWAVE PULSE 

In this Appendix we evaluate analytically and numer- 
ically the initial density perturbations in the ranges of 
laser parameters (amplitude, duration, and the beat fre- 
quency) relevant to the numerical examples of this paper. 
We consider the excitation of EPW at the plasma border 
z = and take the Gaussian temporal profile of laser 1241 
Then, the normalized amplitude e(t) = 6n(t)/no of near- 
resonantly driven electron density perturbation obeys the 
weakly nonlinear equation 

[n~ l d/dt - i5uo(t)]e = -i{A/2) 2 e- 2t2 / T Z , (Al) 

where 5u(t) = 5Cj + (3/16)|e(i)| 2 + (A/2) 2 e - 2t2 ^ , 5u> = 
dui/fl, \5w(t)\ <C 1. We will solve this equation numer- 
ically with the initial condition e(— 00) = (plasma is 
quiescent before the laser pulse arrival). 
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FIG. 9: Temporal evolution of the EPW driven by the two- 
frequency laser pulse of duration ui p tl = 300. The laser inten- 
sity profile |ao(t)| 2 + |ai(t)| 2 is shown with a dotted line. Nu- 
merical solution to Eq. 1A11 gives the EPW amplitude in the 
resonant (Q = uo p , solid line) and near-resonant (Q = 0.9u> p , 
dashed line) cases. 

Note, that Eq. (|A"TJ) can be solved analytically by 
means of the perturbational approach when \8u>i\ > A = 
(l/2)^/3(A/2) 4 (Ref. [13), and 5ut l ^> 1. Then, e(t) 
adiabatically follows the temporal profile of ponderomo- 
tive force, and, in the zero-order approximation, the time 
derivative can be omitted in Eq. (|A1I) . First-order ap- 
proximation is obtained via linearization of Eq. QAlj l with 
\e(t)\ 2 ^ [{A/2) A /5u 2 ]e- At2 / T t. The solution reads 

e{t) « -i(a/2) 2 e iSu;t - inn ^ 

x f e -2T 2 /^-«w+^T L *(r) dTj (A2) 
J — 00 

where $(x) = Kierfc(2x/rx) + K%exic{^/2x /tl), erfc(y) = 
(2/^)Jy C °e- t2 dt, ki = (30F/64)(yl/2) 4 /to 2 , and 

k 2 = ^/8{A/2) 2 . 

Figure illustrates the evolution of the EPW am- 
plitude for the resonant (f2 = lo p ) and near-resonant 
(Q. = 0.9ujp) excitation. In the near-resonant case, 
the beatwave pulse amplitude varies slowly: |<5w7l| « 
30. Resonant plasma beatwave [numerical solution of 
Eq. IjAip ] exhibits oscillations due to the periodic dephas- 
ing produced by the relativistic frequency shift [8|, [10( . 
This effect is eliminated by introducing the detuning 
SCj Ki — 3 A = —0.1 (compare dashed and solid curves 
in Fig. |5J). Under parameters of Fig. [51 Suj dominates 
the relativistic frequency shifts [(3/16)|e(£)| 2 < 0.002 and 
l°o(*)/2| 2 < 0.01], and the numerical solution of Eq . (|A1|) 
almost coincides with the analytic approximation l|A2|l . 

The amplitude of the seed plasma wave from which the 
SFRS can grow inside plasma can be estimated using the 
the amplitude of wake plasma wave left behind the laser, 
\e(t — > +00) I . We evaluate it numerically from Eq. (|A1|) 
and plot as a function of tl for various negative detunings 
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5u> in Figs. ^] (with the maximal laser intensity as of 
Fig. |3) and Fig. ^2 (with the maximal laser intensity as of 
Fig.EJ. For the large negative detunings (5u> < — 0.055w p 
for Fig.HUI and Sco < -0.0125w p for Fig. [TTJ the plasma 
wake amplitude drops with tl monotonously, while for 
smaller magnitudes of \6u>\ (numerical examples are given 
for 8uj = — 0.04w p in Fig. an d for 5uj = — 0.01w p 
in Fig. Ill|) it exhibits an oscillatory behavior. As the 
pulse duration tl grows, the wake amplitude periodically 
reaches the value \e(t — > +oo)| « 0.7 (Fig.EJ and \e(t — > 
+oo)| w 0.5 fFig. lTTJl . This is the signature of relativistic 
bi-stability of the system 0, 0] . 
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FIG. 10: Amplitude of the plasma wake excited by the de- 
tuned heatwave as a function of laser pulse duration; a 2 — 
0.04. 
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FIG. 11: Amplitude of the plasma wake excited by the de- 
tuned heatwave as a function of laser pulse duration; a 2 — 
0.005. 
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